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ABSTRACT 

In March 2009 the well-studied quasar, PG 0844+349, was discovered with Swift to be in 
an X-ray weak state. A follow-up XMM-Newton observation several weeks later generated a 
good quality spectrum of the source, showing substantial curvature and spectral hardening. In 
combination with archival data at two previous epochs when the source was in a bright state, 
we examine the long-term spectral and timing properties of PG 0844+349 spanning nearly 
ten years and a factor of ten in brightness. Partial covering and blurred reflection models are 
compared to the data at each flux state while attempting to maintain consistency between 
the various epochs. In terms of the blurred reflection model, PG 0844+349 is in a reflection 
dominated state during the 2009 X-ray weak observations, which can be understood in terms 
of light bending. Moreover, the light bending scenario can also account for the short-term (i.e. 
~ 1000 s) spectral variability in the source. Other models cannot be decisively ruled out, but 
we note distinguishing features of the models that can be explored for in higher signal-to-noise 
data from current and future observatories. 

Key words: galaxies: active - galaxies: nuclei - galaxies: individual: PG 0844+349 - X-ray: 
galaxies 



1 INTRODUCTION 

The X-ray spectra of active galactic nuclei (AGNs) can exhibit such 
significant differences from one epoch to the next that it can be- 
come very difficult to extract meaningful constraints from single- 
epoch spectra. Multi-epoch data are particularly useful to constrain 
model components that fluctuate on vastly different time scales, 
thereby removing excessive freedom in some fit parameters (e.g. 
Gallo et al. 2010, 2007a,b; Grupe et al. 2008; Zoghbi et al. 2008) 
and allowing for more insightful discussion. 

One attractive proposition is to catch AGN in a particularly 
low X-ray flux state (Gallo 2006). In principle, the low-flux state 
can provide the best perspective of the AGN environment, for ex- 
ample the nature of the accretion disc and associated relativistic 
features (e.g. Miller 2007) and/or the surrounding ionised plasmas 
(e.g. Longinotti et al. 2008) since the direct continuum (i.e. the 
power law component) is substantially suppressed (e.g. Schartel et 
al. 2010, Ballo et al. 2008, Grupe et al. 2008; Vignali et al. 2008; 
Miniutti et al. 2009). Moreover, examining the spectral changes 
from high- to low-flux states can potentially reveal the origin of 
the variability, whether it be: light-bending (e.g. Miniutti & Fabian 



2004), absorption (e.g. Turner & Miller 2009), or beaming (e.g. 
Gallo et al. 2010, 2007b). 

Efforts to identify objects in the low X-ray flux state by mon- 
itoring AGNs with Swift and triggering a pointed XMM-Newton 
observation when an object is discovered to be in a low-state are 
proving advantageous (e.g. Grupe et al. 2007, 2008). We have re- 
cently applied this strategy to catch the type-I quasar, PG 0844+349 
(z = 0.064) in an X-ray weak stateQ 

PG 0844+349 is a particularly interesting radio-quiet AGN. 
Sufficiently bright in the optical to be termed a quasar, it possesses 
strong Fe II emission and weak [O III] like a narrow-line Seyfert 1 
(NLS1), though its FWHM(H/3) (2420 kms" 1 , Boroson & Green 
1992) is slightly above the oft-quoted upper limit for a NLS1 of 
2000 km s~ 1 (Osterbrock & Pogge 1 985). Its black hole mass esti- 



An X-ray weak state implies that the AGN X-ray emission is significantly 
diminished with respect to the AGN's UV flux. This is not necessarily the 
same as a X-ray low flux state, which makes no statement about the UV 
emission and typically refers to the brightness of the source relative to pre- 
vious X-ray observations. The X-ray state can be quantified by a ox as we 
describe in Sect. 3.1. In 2009, PG 0844+349 was in a X-ray weak and low- 
flux state. 
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mated from reverberation mapping is 9.24±3.81 x 10 s M Q (Peter- 
son et al. 2004). Several values can be estimated for its Eddington 
ratio based on the literature (e.g. Grupe et al. 2010; Vasudevan & 
Fabian 2009), but the value is apparently always sub-Eddington. 

Early in the XMM-Newton mission, an observation of 
PG 0844+349 garnered attention as it was found in a histor- 
ically high-flux state compared to previous X-ray observations 
(Brinkmann et al. 2003). Pounds et al. (2003) reported absorption 
features in this same data set that they interpreted as high-mass, rel- 
ativistic outflows. A subsequent reanalysis of the data along with 
a deeper observation by Brinkmann et al. (2006) could not con- 
firm the absorption features reported by Pounds et al. and attributed 
them to calibration uncertainties (see also Vaughan & Uttley 2008). 

In March 2009, PG 0844+349 was observed with Swift as part 
of an unrelated project (PI: Vanden Berk). A quick look at the X- 
ray data revealed the quasar had entered a low-flux state. A Swift 
Target of Opportunity observation (ToO) was initiated in April to 
determine if the AGN was still X-ray dim. As it was found to be so, 
an XMM-Newton ToO was triggered a few days later and it is these 
new data, in combination with the previous XMM-Newton observa- 
tions, that we report on here. 

The observations and data reduction are described in the fol- 
lowing section. In Section 3, the general properties of the data at 
each epoch are explored. In Section 4 we fit the average spectrum at 
each epoch with various physical models to derive a self-consistent 
interpretation for the long-term (yearly scale) changes. We compare 
the model, consider the short-term (intra-observation) variability, 
and discuss our results in Section 5. We summarise our results in 
Section 6. 



2 OBSERVATIONS AND DATA REDUCTION 

Snap-shot observations of PG 0844+349 were conducted with Swift 
(Gehrels et al. 2004) on two occasions in 2009, the details of which 
can be found in Table Q] The X-ray Telescope (XRT; Burrows et 
al. 2005) operated in photon counting mode (Hill et al. 2004) and 
the data were reduced with the task xrtpipeline vO . 12 . 1. 
Source photons were extracted from a circular region of 47" radius 
(corresponding to an encircled energy fraction of 90 per cent) using 
xs elect. The background events were selected from a nearby 
source-free region with radius 236". Net source counts were low at 
both epochs (126 and 63) and a spectrum was only created for the 
higher count spectrum (2009 March). The UV-Optical Telescope 
(UVOT; Roming et al. 2005) observed PG 0844+349 through the 
various filters at both epochs. The UVOT data were reduced and 
analyzed as described in Poole et al. (2008). All magnitudes are 
corrected for Galactic reddening (Eb-v = 0.037; Schlegel et al. 
1998) and are listed in Tabled 

PG 0844+349 was observed with XMM-Newton (Jansen et al. 
2001) in 2000 and 2009. In 2001 the AGN was in the field-of-view 
during the observation of cluster of galaxies Vik 59. A summary of 
all X-ray observations is provided in Table Q] During the 2000 and 
2009 observations the EPIC pn (Striider et al. 2001) was operated 
in full-frame and large- window mode, respectively. In 2001 the pn 
camera did not collect data. The MOS (MOS1 and MOS2; Turner 
et al. 2001) cameras were operated in small-window mode in 2000 
and 2009, and full-frame mode in 2001. The medium filter was in 
place for the first observation and the thin filter for the other obser- 
vations. The Reflection Grating Spectrometers (RGS1 and RGS2; 
den Herder et al. 2001) also collected data during these observa- 
tions, but will not be discussed here. The data were of lower sta- 



tistical quality when PG 0844+349 was in a low-flux state in 2009, 
and in 2001 the source was not detected as it was observed off- 
axis. The 2000 RGS data were reported by Pounds et al. (2003). 
The Optical Monitor (OM; Mason et al. 2001) operated in imaging 
mode and collected data in the U filter in 2000 and various filters 
in 2009 (Table |2j. No OM data were collected in 2001 as the AGN 
was outside the OM field-of-view. 

The XMM-Newton Observation Data Files (ODFs) from all 
observations were processed to produce calibrated event lists using 
the XMM-Newton Science Analysis System (SAS vlO . . 0). Un- 
wanted hot, dead, or flickering pixels were removed as were events 
due to electronic noise. Event energies were corrected for charge- 
transfer inefficiencies. Light curves were extracted from these event 
lists to search for periods of high background flaring. Significant 
background flares were detected in 2000 and 2001, and those pe- 
riods have been neglected. The total amount of good exposure is 
listed in Table Q] Source photons were extracted from an annulus 
around the source to mitigate the effects of pile-up in the central 
region during 2000 and 2001. Pile-up was negligible in the 2009 
observations and a circular extraction region that included all of 
the source was used. The background photons were extracted from 
an off-source region on the same CCD. Single and double events 
were selected for the pn detector, and single-quadruple events were 
selected for the MOS. EPIC response matrices were generated us- 
ing the SAS tasks ARFGEN and RMFGEN. The MOS and pn data at 
each epoch were compared for consistency and determined to be in 
agreement within known uncertainties (Guainazzi et al. 2010). The 
MOS1 and MOS2 spectra were combined to produce a single MOS 
spectrum at each epoch. 

The XMM-Newton X-ray spectra were grouped such that each 
bin had a signal-to-noise ratio of at least 7 in 2000 and 2001, and 
5 in 2009. Spectral fitting was performed using XSPEC vl2 . 5 . 
(Arnaud 1996). All parameters are reported in the rest frame of the 
source unless specified otherwise. The quoted errors on the model 
parameters correspond to a 90% confidence level for one interesting 
parameter (i.e. a A\ 2 = 2.7 criterion). A value for the Galactic col- 
umn density toward PG 0844+349 of 3.39 x 10 20 cm" 2 (Elvis et al. 
1989) is adopted in all of the spectral fits. K-corrected luminosities 
are calculated using a Hubble constant of Hq=70 km s _1 Mpc -1 
and a standard flat cosmology with Qm = 0.3 and 57a = 0.7. 



3 A FIRST-LOOK AT THE DATA 

3.1 The UV/X-ray spectral energy distribution 

The 2000 observation of PG 0844+349 represents a relatively typ- 
ical state for the AGN. The optical-to-X-ray spectral slope (i.e. 
a. ox ; Tananbaum et al. 1979) was approximately —1.4 and was 
within expected values given its 2500 A luminosity (Just et al. 
2007), where we have estimated the 2500 A values directly from 
the shape of the spectral energy distribution (SED) (Grupe et al. 
2010). In 2001 PG 0844+349 was slight brighter in the X-rays, 
but no OM data are available. Assuming the same optical flux as 
in 2000, the 2001 slope would be a ox w -1.3 In March 2009, 
Swift detected a substantially lower X-ray flux, while no significant 
change was measured in the optical. Consequently the measured 
a ox was steeper (a ox ~ —1.86) and the object was in an X-ray 
weak state. In Figure Q] the UV/X-ray SED of PG 0844+349 in 
2000 and March 2009 (i.e. 2009a) is shown. The follow-up Swift 
observation in April 2009 revealed similar fluxes indicating a con- 
tinued X-ray low state so an XMM-Newton Target of Opportunity 
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Table 1. PG 0844+349 X-ray observation log. The observatory and instrument used to obtain the data is given in column 2. MOS refers to the combine MOS1 
and MOS2. The observation ID is provided in column 3. The UT start date and the good-time interval (after removal of background flaring intervals) are 
given in column 4 and 5, respectively. Column 6 provides the approximate total number of source counts in the 0.5 — 10 keV band. The MOS information in 
columns 5 and 6 is for the combined MOS1 and MOS2. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Epoch 


Telescope and 


Observation ID 


Start Date 


Exposure 


Counts 




Instrument 




(year.mm.dd) 


(s) 




2000 


XMM-Newton pn 


0103660201 


2000.11.05 


11940 


23241 




XMM-Newton MOS 






37930 


48777 


2001 


XMM-Newton MOS 


0107860501 


2001.10.08 


127600 


133484 


2009a 


Swift XRT 




2009.03.14 


3618 


126 


2009b 


Swift XRT 




2009.04.26 


1614 


63 


2009 


XMM-Newton pn 


0554710101 


2009.05.03 


11700 


3716 




XMM-Newton MOS 






28870 


2608 



Table 2. The optical/UV observations of PG 0844+349. The observatory and instrument used to obtain the data is given in column 2. The magnitude and 
associated statistical error is provide for each filter in columns 3-8. While the OM and UVOT filter sets are not identical they are comparable (e.g. Grupe et 
al. 2008). The effective wavelength for each filter is given in A and shown in brackets with the UVOT value preceding the OM. All magnitudes are corrected 
for Galactic reddening of i?e-v = 0.037. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


Epoch 


Telescope 


V 


B 


U 


UVW1 


UVM2 


UVW2 






(5468/5235) 


(4392/4050) 


(3465/3275) 


(2600/2675) 


(2246/2205) 


(1928/1984) 


2000 


XMM-Newton OM 






13.92 ±0.01 








2009a 


Swift UVOT 


14.45 ±0.01 


14.79 ±0.01 


13.65 ±0.01 


13.66 ±0.01 


13.61 ±0.01 


13.64 ±0.01 


2009b 


Swift UVOT 


14.49 ±0.02 


14.79 ±0.01 


13.68 ±0.01 


13.67 ±0.01 


13.55 ±0.02 


13.63 ±0.01 


2009 


XMM-Newton OM 


14.54 ±0.01 


14.98 ±0.01 


13.95 ±0.01 


13.71 ±0.01 


13.79 ±0.01 


13.81 ±0.02 



is the difference between the measured and expected a ox ; Just et 
al. 2007). Following Grupe et al. (2010), the Eddington ratio is 
L/LEdd ~ 0.1 during both X-ray weak and normal states. 



3.2 X-ray variability 

Presented in Figure|2]is the 2-10 keV flux of PG 0844+349 from 
all observations over the last 25 years. The light curve shows sub- 
stantial variability and demonstrates that departures into the X-ray 
low-flux state are common for PG 0844+349. Not included in the 
figure are the ROSAT date, which also show a factor of 6 — 10 vari- 
ability at low energies (0.1 — 2.4 keV) in observations separated 
by about 6-months (Rachen et al. 1996; Wang et al. 2000). Similar 
fluctuations are seen in the XMM-Newton data at energies between 
0.5-2keV. 

The X-ray spectral variability from X-ray high- to weak- 
state are depicted in Figure [3] The high-flux states show a rather 
smooth spectrum across the entire energy band. In contrast the 
low-flux state exhibits substantial hardening and at lower ener- 
gies PG 0844+349 is an order of magnitude dimmer than in the 
high-flux state. The difference spectrum between the high (2000) 
and low (2009) states is computed revealing the shape of the vari- 
able component. This spectrum is fit well with a power law (r = 
2.4±0.03) between 1 — 10 keV, but deviates significantly when the 
fit is extrapolated to lower energies (Figure[4]l. The variable compo- 
nent cannot be fit by a single power law across the entire band, thus 
indicating there is at least another component contributing to the 
long-term variability. In combination with the ROSAT analysis by 




Figure 1. The UV-to-X-ray SED of PG 0844+349 at two separate epochs. 
In 2000 (black circles) XMM-Newton observed the AGN in a normal flux 
state. In March 2009 (2009a, red triangles) Swift observed PG 0844+349 in 
an X-ray weak state. The 2009a SED is representative of the spectral shape 
as it was observed again in April and May 2009. There was only one optical 
filter (U) used in 2000. 



(ToO) was triggered. XMM-Newton also caught PG 0844+349 in 
an X-ray weak state. The measured UV-to-X-ray slope during this 
observations was a ox ~ —1.80 and Aa ox = —0.39 (where Aa ox 
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Figure 2. The 2-10 keV light curve of PG 0844+349 spanning nearly 
25 years. The instruments used to make the measurements are indicated 
beside the data points. The XMM-Newton and Swift fluxes come from this 
work and the Chandra values come from Shu et al. (2010). The earlier flux 
values are taken from Wang et al. (2000). Departures into a low X-ray flux 
state are not uncommon for PG 0844+349. Error bars are included for all 
points except the for the Chandra data where none were reported. The time 
axis is continuous, but scales differently in each panel in order to enhance 
clarity. 
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Figure 3. The spectra from the three XMM-Newton observations are shown 
after correcting for the effective area of the detectors. The MOS and pn 
data are shown with the same colour at each epoch to minimise confusion. 
PG 0844+349 is in a low-flux state in the 2009 and exhibits significant cur- 
vature in its spectrum (blue). The different epochs are identified by colour: 
2000 (black); 2001 (green); and 2009 (blue). 



Wang et al. (2000), the low-energy flux variations may be generally 
more significant in PG 0844+349 than at high energies. 

The combined EPIC light curve from each XMM-Newton ob- 
servations is shown in Figure[5] Flux variations on rapid time scales 
(e.g. < 1000 s) are present in the longer and higher signal-to-noise 
light curves composing the two high state observations. During the 
low-flux state the AGN is constant on the ~ ±15 per cent level. 

Spectral variability is examined for in the longest observation 



Figure 4. The difference spectrum between the high-flux (2000) and low- 
flux (2009) data is fit with a power law (corrected for Galactic absorption) 
between 1 — lOkeV. The fit is then extrapolated over the 0.5 — lOkeV 
band. 
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Figure 5. The 0.2 — 10 keV EPIC light curves from the three observations 
of PG 0844+349 in 100 s bins. The counts from all EPIC instruments that 
collected data at each epoch have been added together. N.B. Not all three 
cameras collected data at each observation so the count rates on the vertical 
axes are not directly comparable between epochs (see Table[T}. Zero on the 
time axes marks the start of each observation. Note that the duration of each 
observation is different. 



(2001) and in the low-flux state (2009) by calculating the fractional 
variability (F var ) at various energies (e.g. Ponti et al. 2004). In 2001 
PG 0844+349 was in a high-flux state corresponding to the bright- 
est smooth spectrum (green data) in Figure[3] The spectral variabil- 
ity is modest, but there is a trend of increasing amplitude with en- 
ergy. The 2009 data are of lower quality as the variability is small 
and the observation is short. There is no detectable spectral vari- 
ability within about ±15 per cent. The F var spectra are presented 
and discussed further in Section|5]and Figure [Tol 

Hardness ratio variability curves in various energy bands dur- 
ing the 2001 epoch show modest to low fluctuations with time. 
There also appears to be no significant correlation between hard- 
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Table 3. The measured energy and flux of an intrinsically narrow (a = 
1 eV) Gaussian profile that is added to the power law model in the 2.2 — 
lOkeV band. Flux units are erg cm -2 s — 1 . The final row shows the im- 
provement to the power law fit when the two free parameters are added. 
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Figure 6. The residuals (data/model) remaining from fitting each spectrum 
with a power law absorbed by Galactic absorption in the 2.2 — 4.5 and 
7 — 10 keV band, and then extrapolating over 0.5 — 10 keV. An excess at 
low-energies (< 2keV) is prominent in all flux states (i.e. the soft excess). 
An excess at high-energies (4 — 8 keV) is more prominent in the low-flux 
state. The inset is an enlargement of the 4 — 8 keV region and the axes are 
the same as on the large plot. The vertical dashed line marks the rest-frame 
6.4 keV. Colours are as described in Figure|3] 



ness ratio fluctuations and count rate. Similar conclusions were 
reached by Brinkmann et al. (2006). 



3.3 Phenomological X-ray spectral model 

In this section we attempt to fit the spectra of PG 0844+349 with 
phenomological AGN models. Fitting a power law modified by 
Galactic absorption to the 2.2 — 4.5 keV and 7 — lOkeV band, 
and then extrapolating over the 0.5 — lOkeV band generates the 
residuals seen in Figure [6] In the high-flux states the power law 
describes the data well above ~ 2 keV and positive residuals are 
seen above the power law below ~ 2 keV (i.e. the soft excess). In 
2009 the photon index of the power law fit is much flatter then in 
the high-flux state and the soft excess appears much stronger. There 
are also significant excess residuals between 5 — 7 keV owing to the 
prominent spectral curvature in the low-flux state. 

The soft excess is well fitted with a black body in addition to 
the underlying power law. The best-fit temperature is comparable in 
the high-flux states: 110±5 eV in 2000 and 107±3 eV in 2001; but 
significantly higher in 2009 (133 ± 6 eV). All the temperatures fall 
in the typical range exhibited by unabsorbed AGN (e.g. Gierlinski 
& Done 2004; Crummy et al. 2006). The measured temperature 
is not correlated in an obvious manner with the black body lumi- 
nosity (e.g. the highest temperature corresponds to the intermediate 
luminosity) and a standard thermal accretion disc origin for the soft 
excess is unlikely. 

A narrow (a — 1 eV) Gaussian profile is added to the fit 
to search for Fe Ka emission at ~ 6.4 keV arising from distant 
material (e.g. the torus). A marginal improvement over the power 
law model is seen at each epoch including such a feature, and the 
measured energies and fluxes are consistent at all epochs. Fitting a 
common feature at all epochs provides the most significant detec- 
tion (Table [3}. When permitting the width of the line to vary, the 
fit prefers a broad feature likely compensating for spectral curva- 



(1) 

Parameter 


(2) 
2000 


(3) 
2001 


(4) 
2009 


(5) 
Combined 


E„ st (keV) 
logF 
Ax 2 


6.34 ±0.10 
-13 41+ - 20 
7.7 


6.44 ± 0.06 
-13 41 +0 - 20 

±0 -** -0.38 

7.3 


6 41+ - 20 
9.5 


6.41 ± 0.06 

-13J58±S;» 
20.7 



ture and/or emission due to He- and H-like iron. While these high- 
ionisation lines may very well be present, they would not account 
for the excess flux at E < 6.4 keV. Our measurements of energy 
and flux are consistent with the recent analysis of Shu et al. (2010) 
who detected a significant 6.4 keV emission line in one-third of a 
Chandra HETG observation of PG 0844+349. 



4 MULTI-EPOCH SPECTRAL ANALYSIS OF 
PG 0844+349 

The spectral changes from high- to low-flux state in PG 0844+349 
are significant (e.g. Figure |3j- In this section we will attempt 
to model the entire data set with more physical models and 
in a self consistent manner. In shape, the low-flux spectrum of 
PG 0844+349 is reminiscent of other AGNs caught in a low-flux 
state (e.g. 1H 0707-495 and Mrk 335). The spectrum of such ob- 
jects being modelled by blurred reflection or some form of partial 
covering. All models include two narrow Gaussian profiles with 
fixed energies and widths (a — 1 eV) that account for distant 
Fe Ka (E — 6.4 keV) detected here and with Chandra (Shu et al. 
2010), and the undetected but expected Fe K/3 (7.058 keV) emis- 
sion line. The normalisation of the Fe K/3 feature is fixed to 0.17 
that of the Fe Ka. All models and parameters are shown in Table|4] 
Based on the statistically best-fit model (i.e. the ionised disc re- 
flection) the 0.5 - 2 keV luminosities of PG 0844+349 (corrected 
for Galactic column density) are ~ 5 x 10 42 and 10 44 ergs _1 
in the low- and high-flux state, respectively. The corresponding 
2 — 10 keV luminosities in the low and high state are w 10 43 and 



7 x 10 43 erg s 1 , respectively. 



4.1 Ionised disc reflection 

Reflection from an ionised disc blurred for relativistic effects close 
to the black hole is often adopted to describe the origin of the soft 
excess (Ross & Fabian 2005; Ballantyne et al. 2001), and has been 
successfully fitted to the spectra of unabsorbed AGN (e.g. Fabian 
et al. 2004; Crummy et al. 2006). In particular, the spectral com- 
plexity seen in the X-ray weak state can be attributed to a reflection 
dominated spectrum. That is when the direct power law component 
seen by the observer is significantly diminished due to obscuration 
(e.g. Fabian et al. 2002) or light bending (Miniutti &Fabian 2004). 

An intrinsic power law continuum plus blurred (kdblur) re- 
flection model is fit to the data of PG 0844+349. The use of multi- 
epoch data permits linking of parameters that are not expected to 
vary over observable periods, for example: disc inclination (i), iron 
abundance (Af £ ), outer disc radius (R ou t)- To further simplify 
the fitting process the inner (Ri n ) and outer disc radius are fixed 
at 1.235 r g and 300 r g , respectively (lr g = GM/c 2 ). Allowing 
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Table 4. The best-fit model parameters for the multi-epoch spectral modelling of PG 0844+349. The model, model components and model parameters are 
listed in Columns 1, 2 and 3, respectively. Each subsequent column refers to a specific epoch. The inner (Ri„) and outer (R ou t) disc radius are given in units 
of gravitational radii (1 r g = GM/c 2 ). The absorber covering fraction (Cf) is given in percentage. Values that are linked between epochs appear in only one 
column. The superscript / identifies parameters that are fixed. Fluxes are corrected for Galactic absorption and are in units of ergcm~ 2 s _1 . 



(1) 

Model 


(2) 

Model Component 


(3) 

Model Parameter 


(4) 
2000 


(5) 
2001 


(6) 
2009 


Blurred reflection 


Power law 


r 

log F 5 _ 10keV 


2.31 ± 0.03 
-11.084 ±0.001 


2.24 ± 0.02 
-10.933 ±0.001 
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Figure 7. Upper panel: The ionised reflection model used to fit the high- 
state (2001) and low-state (2009) spectra of PG 0844+349. The 2000 model 
is very similar to the 2001 and is not shown for clarity. The dotted curve 
and dashed curve are the power law and reflection components, respec- 
tively. The solid curve is the combined model. The reflection component 
dominates the X-ray band in the low-state. Lower panel: Spectral residu- 
als (in sigma) based on this model (see text and Table|4]for details) for the 
high-state (top) and low-state (bottom). 



these parameters to vary reproduced similar values, but did not sig- 
nificantly improve the fit quality. The only blurring parameter al- 
lowed to vary at each epoch is the emissivity profile of the disc, 
which is a power law in radius (J(r) oc r~ a ). The ionisation pa- 
rameter of the disc (£ = L x /nr 2 , where n is the hydrogen number 
density and r is the distance between the ionising source and the 
material), power law continuum (r and normalisation) and normal- 
isation of the reflector are permitted to vary at each epoch. 

The model provides a good fit to the spectra at all three epochs 
(X?/dof = 1.05/474; Figure |7J, and can be considered self- 
consistent. While there are several parameters that fluctuate from 
epoch-to-epoch the primary driver is the diminishing of the power 
law normalisation from high- to low-flux. Indeed, during the X- 
ray weak state the X-ray spectrum of PG 0844+349 is dominated 
by the reflection component (Figure [7J upper panel). We estimate 
the reflection fraction (R) at each epoch by comparing the ex- 



Figure 8. Upper panel: Spectral residuals (in sigma) from modeling the 
high-state (top) and low-state (bottom) of PG 0844+349 with a neutral dou- 
ble partial covering (see text and Table|4]for details). Lower panel: Same as 
upper panel, but for the ionised double partial covering model. 



pected model above 10 keV (where ionisation and absorption ef- 
fects are minimum) with the pexrav model in XSPEC (Magdziarz 
& Zdziarski 1995). During the high-flux state R « 0.8, commen- 
surate with an isotropic emitter. In the X-ray weak state, R » 1 
indicating that extreme effects, such as light bending, may be at 
play. 

The models also show a flattening of the power law slope from 
high- to low-state. While some variability is expected in the shape 
of the power law with time, the change measured here is rather 
significant and difficult to reconcile. Similar behaviour was noted 
in another AGN, Mrk 79 (Gallo et al. 2010). 



4.2 Neutral partial covering 

We consider a scenario where the intrinsic power law continuum is 
partially obscured by a neutral absorber along the line-of-sight (e.g. 
Grupe et al. 2004, Gallo et al. 2004). In this picture, the observed 
spectrum is the combination of direct emission from the power 
law emitter and a highly obscured component. Since a significant 
amount of time passes between observations, we allow for both in- 
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Figure 9. The neutral partial covering model assuming a constant primary 
power law (red dashed line) at all epochs. The changes in the observed 
spectrum are due to variations in the absorber covering fraction and column 
density. For clarity, only the 2000 and 2009 models are shown. 

trinsic continuum changes (i.e. changes in the power law shape and 
luminosity) and changes in the characteristics of the absorber [i.e. 
column density (Nu) and covering fraction (C/)]. The simplest 
case of a single absorber results in a poor fit(x2/dof = 2.09/479) 
and supports the addition of a second partial covering absorber (i.e. 
neutral double partial covering). The double partial covering model 
is akin to two distinct absorbing regions or a column density gradi- 
ent along the line-of-sight (e.g. Tanaka et al. 2004). 

The fit is significantly improved (xS/dof = 1.16/473) with 
the addition of the second absorber. The model parameters and fit 
quality are shown in Tableland Figure[8](top panel), respectively. 
The two absorbers have similar characteristics during the high-state 
observations (2000 and 2001) and together obscure about 70 per 
cent of the power law light. To account for the low-flux state the 
column density and covering fraction of both absorbers increase, 
and together diminish the source flux in the 0.5 — 10 keV band by 
about 93 per cent. We also note that the intrinsic power law flux 
dropped by about 67 per cent in 2009 from its 2001 high. Note- 
worthy is the particularly steep spectrum that is required for the 
continuum (r ~ 3) at each epoch. 

A statistically worse, but not unreasonable fit, can be obtain 
by maintaining the shape and luminosity of the primary power law 
constant at all three epochs (xt/doi — 1.26/477; Figure[9}. In 
this case, variations in the nature (covering fraction and column 
density) of the absorber account for the long-term variations in 
the shape and flux of the spectrum. The intrinsic luminosity of the 
power law in the 0.5 — lOkeV band is 5.2 x 10 44 ergs^ 1 . The 
power law photon index remains steep (r = 2.93 ± 0.03). 

4.3 Ionised partial covering 

A modification of the partial covering model discussed above is 
to consider an absorber that is ionised (Reeves et al. 2008) rather 
than neutral. The ionised material will preferentially absorb inter- 
mediate energy X-rays giving rise to the different spectral slopes 
on either side of E ~ 2 keV. 

A single ionised absorber partially covering the primary con- 
tinuum is a considerably better fit than the single neutral absorber 



(X^/dof = 1.36/476), but significantly worse than the neutral 
double partial covering and ionised reflection models. The addi- 
tion of a second ionised absorber was attempted to improve the fit. 
In initial fits the ionisation parameter of the second absorber was 
comparable at all three epochs so it was linked in subsequent fits. 
The improvement over the single ionised absorber was substantial 
(Ax 2 = 160 for 7 additional parameters). The fit residuals are dis- 
played in Figure[8] Notably, the ionisation of the absorber does not 
follow the predicted flux of the continuum as would be expected if 
the power law is the sole source of ionisation for the absorber and 
if the absorber was of high density. 

The two absorbers collectively obscure about 50 per cent of 
the power law continuum in the high-state and about 90 per cent in 
the low-flux state. The moderately low ionisation parameters pre- 
dicted with this model will produce a spectrum rich in absorption 
features that should be detectable in high-resolution spectra even in 
the high-state. An analysis of the 2000 XMM-Newton data (includ- 
ing the RGS) by Pounds et al. (2003) resulted in claims of highly 
blue-shifted absorption lines due to a high- velocity outflow, but was 
later dismissed by Brinkmann et al. (2006). 

Similar to the neutral partial covering, the long-term variations 
could be modelled with a constant power law and varying absorp- 
tion parameters. The fit is slightly worse than that presented above 
(X^/dof = 1.16/473) and the intrinsic power law is again steep 
(r = 2.77 ±0.30). 



5 DISCUSSION 

In the X-ray weak state PG 0844+349 shows considerable spectral 
hardening and curvature compared to the relatively smooth spec- 
trum during the bright state. Blurred ionised reflection models fit 
the multi-epoch spectra very well in a self-consistent manner. Par- 
tial covering models, with either neutral or ionised absorbers, also 
work reasonably well. The reflection and absorption models appear 
spectroscopically similar in the 0.5 — 10 keV band, but each model 
appears different at higher energies and predict different temporal 
behaviour. 

In the blurred reflection model the principle difference be- 
tween the high- and low-state is the prominence of the direct power 
law continuum reaching the observer. There is more than a factor 
of 10 difference in the 0.5 — lOkeV power law flux between the 
high- and low-state. In fact, during the low-state the reflection com- 
ponent dominates the EPIC spectrum and the measured reflection 
fraction is R >> 1, The attenuation of the direct power law flux 
could be attributed to obscuring of the component (e.g. Fabian et 
al. 2002) or by light bending effects (Miniutti & Fabian 2004). 

In addition, the reflection models also indicate variations in 
the flux and ionisation of the reflection component. These changes 
are manifested in a model-independent way in the difference spec- 
trum (Figure 0- Rather than being a simple power law (character- 
istic of normalisation changes alone), the difference spectrum has a 
soft excess. Ionisations changes in the reflection component would 
modify the soft excess. 

The light bending scenario makes specific predictions about 
the rapid temporal behaviour of the power law and reflection com- 
ponents in various flux states that can in turn be tested. In exam- 
ining the 2001 F var spectrum (Figure [Tol left panel) we consider a 
situation where the power law component fluctuates in brightness 
and dominates the variability (i.e. the reflection component varies 
very little in comparison). Such conditions are expected from the 
light bending model when the AGN is in the bright state (Mini- 



PG 0844+349 in an X-ray weak state 9 




0.5 1 

Observed energy (keV) 




0.5 1 

Observed energy (keV) 



10 



Figure 10. The normalised rms spectrum calculated for the 2001 (left panel) and 2009 observations (right panel). On the assumption that only the power 
law normalisation fluctuates, the expected i*Var based on the blurred reflection (black solid), double neutral partial covering (red dashed) and double ionised 
partial covering (blue dotted) models are shown. The single ionised absorber produces a similar F vaT as the double ionised absorber and is not shown. The 
light curves used are in 1000 s bins. 



utti & Fabian 2004). Such a model agrees well with the data of 
PG 0844+349 (Figure Qojleft panel). 

Likewise, the expected F va _ r based on the partial covering 
models is also shown in Figure [10] The models are adequate given 
the signal-to-noise, but both models clearly miss the data in some 
energy bins. The absorption models likely require variability of 
an additional parameter to reconcile the F va , T . Gierlihski & Done 
(2006) consider changes in the ionisation parameter of the absorber 
(in addition to power law flux) to reconstruct the F val of some 
sources that show the recognised peak at intermediate energies. 

In a similar fashion we conducted the same exercise for 
PG 0844+349 in the low-flux state (right panel Figure [lOj. The 
spectrum is of lower quality than the 2001 observation as the am- 
plitude of the variations is low and the observation is short, and 
consequently no model can be ruled out with statistical certainty. 
However, we note that the models predict significantly different be- 
haviour and with a longer low-flux observation (and higher signal- 
to-noise) the models could potentially be distinguished. 

If due to partial covering, the simultaneous UV/X-ray data 
could potentially place some restrictions on the geometry of the 
system depending on the nature of the absorber. The UV observa- 
tions (all longward of the Lyman limit) do not appear to exhibit 
diminishing flux as the X-rays do. The fact that we do not see UV 
variability simultaneous with the X-rays implies that any partial 
covering absorber, if present, has to be dust-free and confined to 
the X-ray emitting region. This requires the existence of very dense, 
cold blobs that are a substantial fraction of the source size, occupy- 
ing a small region close to the source. Ostensibly this would imply 
the presence of a strong iron fluorescence line that is not observed 
in the data. Such a strong feature would be detectable unless the 
covering fraction was small or in a rather contrived geometry (see 
Reynolds et al. 2009 and Miller et al. 2010 for opposing opinions). 

The UV spectrum of PG 0844+349 also displays relatively 
weak Civ absorption (Brandt et al. 2000). The feature is expected 
to be much stronger if the X-ray weak state of PG 0844+349 were 
due to absorption from a BAL-related phenomenon (Brandt et al. 



2000). UV spectroscopy of PG 0844+349 when it is in a X-ray 
weak state would confirm the absence of a BAL absorber. 

The inferred intrinsic photon index measured in PG 0844+349 
is steeper than the canonical value of V ~ 1.9 that is normally 
adopted (e.g. Nandra & Pounds 1994). This seems to be indepen- 
dent of the model used and the flux-state of the AGN (except for the 
2009 blurred reflection interpretation, see Sect. 14. It . There are po- 
tential explanations for this, for example if due to Comptonisation 
PG 0844+349 could have a cooler corona than the average AGN 
leading to a steeper spectrum. It is difficult to put this in to context 
of the NLS1/BLS1 nature of PG 0844+349. NLSls are known to 
exhibit steeper spectra than BLS1 (e.g. Boiler et al. 1996; Brandt 
et al. 1997), but this is only based on the observed spectrum that is 
modified by reflection and warm absorption. A true comparison of 
the intrinsic NLS 1 and BLS 1 photon index to examine for physical 
difference between the two classes is still lacking. However, in the 
near future, as BAT detections become more significant and hard 
X-ray imaging becomes possible, we will be in position to answer 
this question. 

Each model has unique spectral signatures that could be dis- 
cerned with higher spectral resolution or broader energy bandpass. 
The ionised absorber predicts absorption features in the spectra at 
6 — 8keV due to various transitions in iron. Such features may 
be lost in the EPIC data due to the moderate spectral resolution. 
However, Shu et al. (2010) present the Chandra HETG spectrum 
of PG 0844+349 when it was in a comparable flux state as the 2000 
and 2001 XMM-Newton observations. The HETG spectrum is of 
modest quality, but shows no indication of absorption features in 
the iron region (see figure 1 of Shu et al. ). A higher signal-to-noise 
HETG spectrum could potential constrain such features. 

Finally, all three models make vastly different predictions of 
the spectral shape and flux above 10 keV. The absorption models 
have steeper intrinsic spectra (r = 2.4 — 3) and a weaker reflection 
component, whereas the blurred reflection model predicts a much 
flatter spectrum. The reflection model also predicts a 10 — 20 keV 
flux that is twice as large in the high-state and 4 — 5 x greater in 
the low-state than the absorption models. Based on our models, the 
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brightest flux reached by PG 0844+349 in the 14 - 195 keV band 
is ~ 1CP 11 erg cm -2 s" 1 and fluctuates to lower fluxes. This peak 
value is still a factor of 2 lower than the current detection limit of 
the Swift BAT survey (Tueller et al. 2010). 

Through multi-epoch X-ray observations the reflection- 
dominated low-state spectrum of PG 0844+349 appears more 
clearly. Similar observations of several other AGN also support 
the notion that the AGN X-ray low-flux state may be due to a di- 
minished power law that in turn reveals the underlying reflection 
spectrum. Multi-epoch observations are necessary to distinguish 
between competing models. Placing constraints on the nature of the 
long-term spectral variability is educative and feasible with current 
instruments. 



6 SUMMARY 

A 15 ks XMM-Newton ToO observation of the Seyfert 1 
PG 0844+349 was triggered when it was discovered with Swift in 
2009 that the AGN was in an X-ray weak state. Here we conducted 
an analysis of the X-ray spectrum and the long-term variability by 
combining data of the AGN from previous epochs when it was in a 
bright state. 

PG 0844+349 is fitted with various physical models (blurred 
reflection and partial covering). While no model is conclusively dis- 
missed the blurred reflection model nicely describes the long- and 
short-term variability in a consistent manner. The two classes of 
models predict distinguishing characteristics that could be revealed 
with deeper observations (i.e. high signal-to-noise and broader en- 
ergy band) of PG 0844+349 in the low-flux state with current X- 
ray observatories. In the future, such observations will be much 
simpler. The high-energy imaging capabilities and calorimeter res- 
olution of Astro-H and 1XO will provide the opportunity to simul- 
taneously measure the broadband spectrum while achieving high 
spectral resolution below 10 keV, and NuSTAR will be the first to 
image the high-energy band. 
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